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Abstract 

We consider a one-dimensional directed polymer in a random potential which 
is characterized by the Gaussian statistics with the finite size local correlations. It 
is shown that the well-known Kardar's solution obtained originally for a directed 
polymer with (^-correlated random potential can be applied for the description of 
the present system only in the high-temperature limit. For the low temperature 
limit we have obtained the new solution which is described by the one-step replica 
symmetry breaking. For the mean square deviation of the directed polymer of the 
linear size L it provides the usual scaling ~ aL^'' with the wandering exponent 
C = 2/3 and the temperature-independent prefactor. 



1 Introduction 



In a wide variety of physical systems one is interested in the behaviour of a fluctuating 
linear object (with finite line tension) interacting with a quenched random potential. 
The object under consideration may be a dislocation in a crystal, a domain wall in a 
two-dimensional magnet, a vortex line in a superconductor, a fiuxon line in an extended 
Josephson junction and so on, but following Ref. this class of problems is traditionally 
discussed in terms of a directed polymer in random media or symply " directed polymer" . 

Quite naturally the best understanding has been achieved for the simplest one- 
dimensional case when the displacements of a directed polymer can occur only in one 
direction. In such case a directed polymer in continuous approximation can be described 
by the Hamiltonian 




where J is the linear tension, t is the longitudal coordinate (0 < t < L) and x{t) is the 
transverse displacement of a polymer with respect to a straight line. The simplest (or 
maybe one should better say the most easily treatable) assumption on the distribution 
of random potential v{x,t) consists in taking it to be Gaussian with 



v(a;,t) = 0; v{x,t)v{x' ,t') = 2V{x - x')5{t - t'). (1.2) 

Here and further on an overbar denotes the average over the realizations of quenched 
random potential. 

Let us assume that at t = the position of a polymer is fixed: x(0) = 0. Then 
the quantity of interest is the typical deviation of the polymer "trajectory" x{t) from 
the origine. More precisely, one would like to know the dependence on L of the average 
square deviation of the polymer at the ending point x = L, which in the limit L — > oo 
is expected to be described by the simple scaling: 

(a;2(L)) ~ aL^^ (1.3) 

where angular brackets denote the average over thermal fluctuations, and ( is the so- 
called wandering exponent. 

In the absence of the random potential the situation is trivial and the wandering 
exponent ( is equal to 1/2. In this case the trajectory deviates from the origine only 
due to the thermal fluctuations, the prefactor a in the scaling law ( |1.3|) is proportional 
to the temperature, so at zero temperature (x^(L)) = 0. 

In the presence of a quenched random potential the situation is getting much more 
complicated. Now besides the thermal fluctuations, the trajectory is pushed away from 
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the origine also due to the randomness in the background potential landscape, so that the 
scaling law (|1.3| ) could be governed by a new non-trivial wandering exponent. Moreover, 
since for a generic random potential the ground state trajectory of the Hamiltonian ( |1.1|) 
typically drifts away from the origine, the scaling law ( |1.3D can be expected to hold also 
in the zero temperature limit. 

It is widely believed that for a whole class of locally correlated random potentials 
(such that the function V{x) in (|1 . 2| ) is quickly decaying for \x\ — *• oo) the wandering 
exponent ( is universal and equal to 2/3. This conclusion is based on the finite tem- 
perature exact results of Refs. @, Q and has been also confirmed by zero temperature 
numerical simulations of the discrete version of the directed polymer problem (see 
also Ref. Q for later references). However, both the calculation based on the reduc- 
tion to damped Burgers' equation with conservative random force and Bethe ansatz 
calculation in terms of replica representation 0] are valid only for the case of stricktly 
(5-functional correlations of random potential. 

In this work we present an attempt to generalize the second of these solutions to 
a more physical situation when correlations of random potential have finite correlation 
radius (in transverse direction). 

In Sec. 2 we review some details of the approach developed by Kardar in Ref. |^ 
for 5-functional correlations of random potential and use it to find the temperature 
dependence of a prefactor in Eq. ( |1.3| ) which turns out to be of the form 

a oc (1.4) 

In Sec. 3 the applicability of Kardar's solution for the approximate description 
of the system with the finite radius r of the random potential correlations is discussed. 
We demonstrate that it can be used only in the high temperature limit T ^ Tq oc r^/^, 
whereas at low temperatures the solution has to have essentially different structure which 
in principle could lead to the change in wandering exponent. The impossibility to apply 
the description with the help of the Kardar's solution at arbitarily low temperatures 
follows already from Eq. ( |1.4| ). One can expect that if in the low temperature limit 
the typical trajectory goes away from the origine its drift should be determined by 
the quenched random potential and not by the effects of the thermal fluctuations. It 
means that in the system with reasonable short-scale regularization the divergence of 
the prefactor suggested by Eq. (|1.4|) at low temperatures can be expected to saturate. 

In Sec. 4 the low temperature solution of the regularized problem is found for the 
particular choice of the random potential correlation function V{x). The form of this 
solution can be described in terms of the effective one-step replica symmetry breaking 
ansatz. In this case one recovers the scaling law ( p..3| ) with the same wandering exponent 
C = 2/3 but with a temperature-independent prefactor. 
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2 The solution of the unregularized problem 



2.1 The relation between exponents 

The idea of the Kardar's approach is based on the indirect calculation of the wan- 
dering exponent ( by analyzing the scaling of the typical sample to sample fluctuations 
of the free energy. Suppose that the typical fluctuations of the free energy (produced by 
the random potential) scale as 

6F oc L'^ (2.1) 

where the exponent u is known. On the other hand, if the typical deviation of the 
trajectory from the origine is equal to x, then the loss of the energy due to the elastic 
term in the Hamiltonian ( |1.1|) must be of the order of Jx^/ L. Ballancing the two energies, 
one can write the following estimate: 

~ ^ oc L-+^ (2.2) 



Then, according to the definition of the wandering exponent ( |1.3| ), one finds the following 
simple relation between the two exponents: 



2C = tu + l (2.3) 



2.2 The rephca method 



The scaling of the free energy fluctuations with the size of the system L can be relatively 
easy investigated in terms of the replica method . To this end one has to calculate 
the average 



Z{n) = (2.4) 
of the n-th power of the partition function: 

ZM = / Z)xWexp|-»i| (2.5) 

Jo<t<L [ I J 

obtained by the integration over all the trajectories with x{0) = 0. 
According to the definition of the free energy: 

F=-T\nZ[v] (2.6) 
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the replica partition function Z{n) can be represented as follows: 



Z{n) = exp 



T 



[2.7) 



The average in Eq. (|2.7| ) as everywhere above is calculated over the realizations of 
random potential v{x, t). On the other hand the free energy F = F[v] is itself the sample- 
dependent random quantity, whose distribution function we shall denote as P{F). Then 
Eq. (p^) can be rewritten as 



Z{n) = J dFP{F) exp 



T 



(2.8) 



which is nothing else but the Laplace transform of the free energy distribution function 
P(F). 

The free energy corresponding to the replica partition function ( p.4|) can be naturally 
defined as: 



F{n) = -T\nZ{n) (2.9) 

Although this quantity can be calculated only for integer n, according to the standart 
ideology of the replica approach it has to be considered as a function of the continuous 
parameter n which implies a necessity of an analytic continuation in n. 

Let us represent the free energy F{n) of the replicated system as a series in powers 
of the replica parameter n: 

oo TP 

k=i ■ 

Then, taking the k-th derivative over n at n = from both sides of the Eq.( |2.9| ) for the 
k-th order of the free energy fluctuations one finds: 

p^^pi-k^ (2.11) 
where double overbar denotes the irreduceble average. 



2.3 The Bethe ansatz type solution 

In the framework of replica approach the statistical meachanics of the irregular system 
is analyzed by considering the statistical mechanics of the regular system in which the 
disorder manifests itself in the form of the interaction between n identical replicas of 
the original system. For the system described by the Hamiltonian (|1 . 1| ) and the Gaus- 
sian statistics of the random potential, Eq. (|1.2| ), the averaging of Z^[v\ over disorder 
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leads to the expression for Z{n) the form of which corresponds to the following replica 
Hamiltonian: 



Hre,i = 1 dt\^T. [-^j - ^ i:VMt) - x,{t)] I . (2.12) 

Eq. ( |2.12| ) has a form of the Eulidean (imaginary time) action describing the quantum- 
mechanical system of n particles with the mass J/T and interaction V{x)/T'^. The same 
system can be described by the quantum-mechanical (operator) Hamiltonian 

rp n -in 

^ = E - ;^ E y^^a - X,) (2.13) 

a=l ^ a,6=l 



which for the classical partition function defined by the Hamiltonion (|2.12| ) plays the 
role of the transfer matrix. 

In the limit of infinite size {L oo) the free energy of a system (for any boundary 
conditions) is dominated by the highest eigenvalue of transfer matrix or, in our case, by 
the lowest eigenvalue Eq of the quantum-mechanical Hamiltonian (|2.13|) 



F{n) = TEo{n)L (2.14) 

For any integer n the lowest eigen- value corresponds to the fully symmetric (nodeless) 
wave-function which for the case of local correlations in x direction 

V{x) = u6{x) (2.15) 
has been found exactly by Kardar [Q: 



'^o[Xa] = exp -33 E - Xb\ (2.16) 

\ a,b=l J 

where for our choice of notation 

ae=^ (2.17) 

The energy of this state is equal to 

^o(n) = - ^n(n^ - 1) (2.18) 

where the first term describes the trivial contribution to E{n) related to the terms with 
a = b in the second sum in Hamiltonian (12. 131). 
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Substitution of Eq. (|2]T8|) into Eq. (^J§ gives 



F(n) = Fin + -F^n^ 

D 



where 



and 



V{0) Ju^ 



T 



6T4 



L 



(2.19) 



(2.20) 



rp4 



L 



Comparison of Eq. (|2.11| ) with Eq. ( p.l9| ) shows that for L 
energy (per unit length) of a random polymer 



(2.21) 



oo the average free 



/ 



lim 

L— >oo 



F{L) 



is given by the linear in n contribution to F{n): 



f 



llim^(^) 



F\ 
L 



V{0) ^ Ju^ 



(2.22) 



(2.23) 



Ln-o n L T ' GT"^ 

in which the first (formally divergent) term always dominates. Therefore the average 
free energy of the system could be defined only after proper short-scale regularization 
of the starting Hamiltonian. 

However, the fluctuations of the free energy are quite well defined without any reg- 
ularization. According to Eqs. ( |2.11| ) and (p.21| ) the typical value of the free energy 
fluctuations can be estimated as: 



/ = X 1/3 / „ X 1/3 ( Ju^y^ 
5F~(|f3|) =(t^|F3|)^ =(^^j L 



Therefore, according to Eq. 
finds the following result: 



'1/3 (2.24) 
for the average square deviation of the trajectory one 



(2.25) 
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3 Introduction of the regularization 



Apparently the divergence of the average free energy, Eg. ( |2.23|) , is removed if one takes 
into account that in a physical system correlations of random potential should be de- 
scribed by a smooth function with a finite correlation radius [and therefore a finite value 
of ^(0)]. However in such case the quantum-mechanical problem defined by Eq. (|2.13|) 
cannot be solved exactly. Nonetheless it seems reasonable to assume that for narrow 
enough V{x) one can still use the expression ( p.23|) in which now u should stand for 



u= / dxV{x) (3.1) 



It is easy to understand that such approximate description [based on Eqs. ( |2.1(j| ) 



( |2.21| )] at low temperatures has to fail. One has to remember that in case of the 6- 



functional interaction the wave function (|2.16| ) is constructed as a generalization of 
two-particle problem wavefunction 

\E'(xi, X2) = exp(-ae|xi - X2I) (3.2) 
On the other hand in the case of a rectangular well: 

^ = 1!^ fhh' (3.3) 

^ ^ [ for |x| > r ^ ^ 

(for which u = 2Vr) the wave function of the two-particle problem for [xi — a;2| > r also 
has the form ( p.2|) with 

-=-jfS) (3.4) 



where 



r \ 



, . j 2:^/2 for z < 1 
'^'^ " . for . » 1 ^3.5) 



and 



To = {JVr^y/'' (3.6) 

It is not hard to check that the condition aer <^ 1 for this wave-function to be wide in 
comparison with the well width coinsides with the condition T ^ Tq. In such case the 
value of se given by Eqs. (p.4|)-( pl6D coinsides with (|2.17|) . In the opposite limit T ^ Tq 
(aer ^ 1) the two-particvle wave-function is almost completely localized inside the well 
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and cannot be used as a building block for the construction of the solution of n-particle 
problem. 

This gives a clear indication that for ar6zirar?/ finite- width form of the function V{x) 
describing the correlations of random potential the applicaton of the Kardar's solution 
for the description of random polymer can work only at high enough temperatures 
whereas in the low temperature limit the solution should be different. 

Quite paradoxically the differentiation of Eq. ( p.23|) shows that for T S> Tq the 
free energy defined by Eq. ( [2. 231 ) corresponds to negative enthropy. One should not be 
too scared of that property since in the approach discussed above the free energy of a 
directed polymer is calibrated in such a way that in the absence of the disorder it is 
equal to zero. Therefore the total free energy will be given by Eq. ( |2.23| ) plus the free 



energy in absence of disorder. This second term will give the positive contribution to 
the enthropy which will overcome the negative contribution from Eq. ( p. 23 ). 



The low temperature solution of the regularized 
problem 



Let us now consider the n-particle problem defined by Hamiltonian ( p.l3|) where 



V{x) = V 



(4.1) 



for < r and is equal to zero elsewhere. Here in comparison with Eq. ( |3.3| ) we have 
introcuded a finite curvature of the potential inside the well which is described by an 
additional free parameter 6 (0 < 6 < 1). At some stage of calculation h will be assumed 
to be much smaller than one. For V{x) of the form (|4.1|) 



u = 2(^-^Vr ^ 2Vr (4.2) 

and therefore the characteristic temperature Tq defining the range of applicability of 
Kardar's approach (T ^ Tq) still can be chosen in the form (|3.6| ). 

The characteristic frequency for small oscillations at the bottom of such truncated 
parabolic well is given by 

2hV 
J2V 

and increases with decrease of T much slower than the depth of the well W = V/T"^. 
Thus the limit of low temperatures may correspond to the case when all particles are 
localized near the bottom of the well. In such limit the ground state energy E[n) for 



^ = \l^2 (4-3) 
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the n-particle system can be rather accurately found by assuming that Eq. (|4.1| ) holds 
for all Xa — Xh. 

In such aproximation the ground state wave-function has a form 



\E' [xa] = exp 
whereas its energy is given by 



1 



AV2n 



a, 6=1 



A 



Q r2 

wJ 



Ein) 



(4.4) 



(4.5) 



(cf. with Refs. and 0). In the following it will be convenient to keep in mind that 
the ratio of Q and W can be expressed as 



(4.6) 



It is not hard to find by a straightforward calculation that for '^{x) of the form (O) 



{{Xa - Xbf) 



A 



1 n 



(4.7) 



'2n ^JImW b 

so for n/W <^ b (that is T <^ 6^/^To) '^(x) is indeed nicely localized at the bottom of 
the well for any integer n and all corrections to E{n) due to non-parabolicity can be 
only exponentially small. 

On the other hand one can easily see why the limit — > is dangerous. The width 
( [4. 71 ) of the wavefunction ( ^.4| ) grows with decrease in n and becomes comparable with 
the width of the well 2r at n ~ 6~^(T/To)^ and therefore for smaller the ground 
state wave-function should have essentially different form. The simplest way to let the 
particles enjoy their mutual attraction while keeping their number in the well not too 
small consists in splittting them into n/k infinitely separated blocks of k particles. The 
energy E{n, k) of such state with broken replica symmetry is given by 



n 



E(n, k) = -Eik) 
k 



n 



-Wk~n- 



k 



and has extremum (maximum) as a function of k. 

Variation of the polymer free energy per unit length: 



(4.8) 



-Wk-Vt- 



f{k)=T 



with respect to k gives an equation for the position of the maximum: 



k 



(4.9) 



10 



the solution of which for Q/W <^ 1 has a form 



fc^~ifli] =i(26)V3jL (4.11) 
2 \WJ 2^ ^ To ^ ^ 



Substitution of Eq. ( [4.11| ) into Eq. ( |4.7| ) then shows that for 



T < fe^/^To (4.12) 
the rephcas belonging to the same block are indeed tightly bound to each other: 

(I)'" (4.13) 

SO the whole picture is really self-consistent. 

Substitution of Eq. (|4.11| ) into Eq. ( [4.9|) gives the temperature independent expres- 
sion: 

f^-lm'^'^ (4.14) 

which shows that in order to find the temperature dependence of / in low temperature 
limit we have to solve Eq. (|4.1(]|) more accurately. That gives 



and 



l/fi\2/3 1/^X^/3 



/^-|(26)V3^ + 1 (26)^3 (4.16) 

/ Jo 4 Iq 

The idea to consider the state in which n replicas are split into n/k infinitely sep- 
arated blocks of k particles has been introduced by Parisi who however applied it 
only to the case of local interaction ( |2.15|) (5-functional correlations in terms of original 
problem) and discovered that free energy as a function of k has extremum at = 0. 
That is equivalent to considering all replicas belonging to the same block right from the 
beginning. Our analysis shows that smearing of the interaction potential leads (at low 
enough temperatures) to the shift of the extremum to non-trivial value of k {0 < k < 1) 
corresponding to splitting of replicas into blocks (that is to replica symmetry breaking). 



Although we have found that in the exremal solution the particles split into k separate 
blocks there are no reasons for these blocks to be infinitely separated from each other. 
The presence of strong attraction between the particles in each block makes it possible to 
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consider such block as a complex particle with the mass kJ/T, the interaction between 
these complex particles being given by k'^V{x)/T'^. The last expression can be expected 
to be very accurate when we consider the temperature interval ( |4.12 ) in which the 



distances between the particles inside each block are much smaller than the well radius 
r. 

Therefore at low temperatures the behaviour of our system in which the particles are 
assumed to be tightly bound mn/k separate blocks can be described by the Hamiltonian 
(12:131) in which 



J ^kJ: 



V{x) A;V(x); u k\ (4.17) 

Our earlier experience tells us that for some values of parameters such system can be 
rather accurately described by the wave function of the form (|2.16| ) in which now Xa 
stand for coordinates of different blocks. The energy of such state will be given by Eq. 
( |2.18|) in which substitutions ( [4.17|) and n ^ n/k have to be made with the first term 
being substituted by Eq. 



E{n, k) 



n 



-Wk-VL 



1 - k 



I'- 



where 




(4.18) 



(4.19) 



is a small parameter if h is small. 

All this leads to appearence in the expression for f{k) of one more term [in compar- 
ison with Eq. (|]9|)]: 



f{k) = T 



-Wk-Vt- 



1 - k 



(4.20) 



which describes the contribution related to mutual interaction between the blocks. Sub- 
stitution of Eq. ( ^.11| ) into the saddle-point equation 



1 _|_ h U/3 



(4.21) 



which is obtained by variation of Eq. ( [4.20|) shows that for B (that is for 6 <^ 1) the 
maximum of f{k) still exists and in the lowest order in h the position of this maximum 
is still given by Eq. ( |4.11 ). 

The applicability of such approach requires that the distances between the blocks 
should be much larger than the size of the well (exactly in same way as when Kardar 
type solution is constructed from the separate particles and not from the blocks): 



ffir ^ 1 



(4.22) 
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Sustitution of Eqs. (|2.17|) , ([4.17|) and ( |4.11|) into Eq. ( [4.22| ) then reduces it to condition 



^ (l - g « 1 (4.23) 

which apparently is equivalent to the same condition 6^1. 

Strictly speaking the expression for f{k) given by Eq. ([4.20|) has also another ex- 
tremum (minimum) at 

k = K,^ (4.24) 



but analogous analysis shows that 



3 / 



ae(A;„)r ^ - - -J ~ 1 (4.25) 

and therfore this second extremum takes place in the domain of parameters where ex- 
pression ( [4.20|) based on the assumption ( [4.22|) can no longer be trusted. 



The form of Eq. (4.18) shows that in the considered case the only non-linear (in n) 



contribution to E{n) is also of the third order in n and corresponds to 

Fs = --^k-L (4.26) 
where we assume 6 -C 1. Using the saddle-point value of the parameter k one gets: 

F,{K) = -{2bf"^L (4.27) 

Correspondingly, for the typical value of the free energy fluctuations one obtaines the 
following temperature independent result: 

5F ~ {T^\F:,\f^ L""" = C^L"^ (4.28) 

where 



Co = (26)2/9(yTo)^/=^ (4.29) 



Finally, for the mean square deviation of the polymer trajectory one again finds (x^) ~ 
aL^/^ with the temperature independent prefactor 
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5 Conclusion 



Thus we have demonstrated that in the case when random potential correlations are 
characterized by finite correlation radius r the solution at low temperatures has es- 
sentially different structure than at high temperatures. Nonetheless the value of the 
wandering exponent in both cases is the same: C = 2/3 . In contrast to the high tem- 
perature limit for which the prefactor a in scaling law ( |1.3| ) is temperature dependent: 
a oc T~^/^, in the low temperature limit it saturates at finite value a oc r^/^. 

Since the value of the wandering exponent for both regimes is the same there are no 
reasons for the sharp transition between these regimes. At very low temperatures the so- 
lution is characterized by the one-step replica symmetry breaking that is the replicas are 
splitted into k well separated blocks. With growth of temperature the distance between 
the blocks becomes comparable with the size of each block. At higher temperatures the 
replica symmetry breaking phenomena can manifest itself only in slight modulation of 
the distance between nearest replicas in comparison with what follows from the "replica 
symmetric" wave- function ([2l6|) . 
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